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ABSTRACT 


The energetics of two numerical models of the ocean 
circulation are studied by use of available potential energy. 
Available potential energy, calculated by an approximation 
equation using a temperature perturbation field, compares 
favorably with that obtained using the exact expression 
uSing a height perturbation field. 

A general circulation model of the ocean is employed to 
study large-scale energy transformations, and a mid-ocean 
dynamics experiment (mesoscale) model is employed to study 
transformations of energy by stationary features and by 
transient motions. Results from the general circulation 
model are consistent with those of Haney (1974). Results 
from the mesoscale model indicate that the transformation 
of energy is from available potential energy to kinetic 


energy of the transient eddy field. 
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IT. INTRODUCTION 


A’“basic part of any study of an atmospheric or oceanic 
circulation system is a study of the behavior of the energy 
in the system and its transformation from one form of energy 
to another. The various forms of energy can be separated 
into three basic forms of energy: kinetic enersy, potential 
energy and internal energy. Haurwitz (1941) showed that 
the potential energy and internal energy, in a vertical 
column, are proportional under hydrostatic equilibrium and 
thus may be combined into a Single parameter, the total 
potential energy. Conversions of energy then occur between 


total potential energy and kinetic energy. 


A. AVAILABLt POTENTIAL ENERGY IN ATMOSPHERIC CIRCULATION 


Lorenz (1955), in a study of atmospheric circulation, 
showed that since total potential energy was three to four 
orders of ASeRuRuEE larger than kinetic energy, transformation 
of energy between the two energy forms was difficult to study. 
in a horizontally stratified atmosphere, heating of the 
atmosphere disturbs the stratification, allowing potential ~ 
energy to be converted to kinetic energy. However, cooling 
of the atmosphere, which removes total potential energy from 
the system, also disturbs the stratification and creates 


conditions for conversion of total potential energy to kinetic 


energy. 
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Lorenz (1955) defined available potential energy as the 
difference between the total potential energy of the whole 
atmosphere and the total potential energy that would exist 
if the mass of the atmosphere were redistributed with 
conservation of potential temperature to form a horizontally 
Stable atmosphere. Then, under adiabatic flow, the sum of 
the available potential energy and kinetic energy is 


conserved. 


B. AVAILABLE POTENTIAL ENERGY IN OCEANIC CIRCULATION 


The use of the concept of available potential energy 
can be carried into the study of oceanic circulations. 
Several studies of the energetics of the ocean circulation 
Systems have been based on available potential energy, as 
in Barrett (1971), and Wright (1972). The study by Wright 
dealt with changes in available potential energy in the 
ocean general circulation. The cooling of cold, high latitude 
water and the surface warming of warm equatorial water 
produces available potential energy by increasing the 
Stratification of the ocean mass. Conversely, the sinking 
of the cold surface water in the north tends to convert 
available potential energy to kinetic energy as it returns 
the water mass to a more stably stratified system. Barrett 
(1971) studied the available potential energy in the transient 
motion of Gulf Stream rings. Results of the study indicated 
that the decay of cold=-core rings in the western Sargasso 


Sea region reduced the available potential energy of the 


1 








region. The presence of decaying cold=-core rings in the 
western Sargasso Sea should have resulted in a rising of 
the mean thermocline in that region. Since the thermocline 
was -not observed to rise significantly in the western 
sargasso Sea region, a local energy source must have been 


present to counteract the energy loss. 


C. AVATLABLE POTENTIAL ENERGY IN NUMERICAL OCEAN 

CIRCULATION MODELS 

This study applied available potential energy calculations 
to two experiments with a numerical model of the ocean general 
circulation; one simulating large-scale, quasi-steady state 
motion and the other simulating mesoscale transient motion. 
It was desired to calculate available potential energy 
using a perturbation height field compared with an approximate . 
poedlacion using a temperature perturbation field. An 
expression for the balance of available potential energy was 
obtained from this approximation equation. The importance 
of the terms in the expression for the balance of available 
potential energy, and their corresponding physical processes, 


could then be examined in the two numerical experiments. 
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fi, DESCRIPRICN OF Toast CAL MODELS 


Piew rere. "experience elle general circulation 
model of the ocean, as described in part A, for large-scale 
quasi-steady state motion. The Mid-Ocean Dynamics Experiment 
model (Robinson, et.al., 1974, MODE, Hot Line News, No. 57.), 
as described in part B, was utilized for the transient motion 


experiment. 


A. GENERAL CIRCULATION MODEL OF THE OCEAN 


The model represents an idealized ocean basin with a 
flat bottom and regular coast. It is everywhere 4 km deep, 
90° longitude wide and extends from 51.25S to 48.75. 

The grid intervals are 3° of longitude and 2.5° of latitude. » 
There are Six vertical levels, placed at 20, 100, 300, 700, 
1500, 3000 meters below the undisturbed surface. The model 
is based on the primitive equations which are integrated in 
time from specified initial conditions. The equations are 
based on the hydrostatic and Boussinesq approximations. 

The ocean is thus assumed incompressible. Salinity is 
neglected and density is a linear function of temperature. 


Tne equations are 


OW 


ges OO. f Uv sing AmV7U ,oru. 
dt Q,acosp oA eS acosP ee alae (1) 
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@2@ Ga- «(7-L)) , 


where the symbols have the following meanings: 


t = time 

A = longitude 

P= latitude 

Z= height 

fi = angular speed of earth's rotation 
f = coriolis force (=2$t sing ) 

a= radius of earth 

@ = acceleration of gravity 

u= eastward velocity component 


v = northward velocity component 


w= vertical velocity component 


T = temperature 
Lae constant reference temperature 
C= density 


14 


(2) 


(3) 


(4) 


(5) 


(6) 





—,. 


= density of water at reference temperature 
= coefficient of thermal expansion 
H = depth of ocean basin 


O{T)= convective adjustment of T if lapse 
rate becomes unstable, 


The lateral boundary conditions on heat and momentum 
at all four walls of the basin are those of perfect insulation 
and zero slip. The ocean floor is thermally insulated, but 
free slip is allowed. Due to the flat bottom, the vertical 
velocity is zero there. The circulation is driven and 
controlled by the wind stress and heat flux at the ocean 
surface. The wind stress is larger in the southern hemisphere 
and meridional wind stress is neglected. The net downward 


heat flux at the ocean surface is expressed as 


@ = on" = Te) } 


where T. is the ocean surface temperature, Ths is an 4) 
atmospheric equilibrium temperature approximately equal to 

the air temperature at 10 meters and Q> is a coupling 
coefficient determined from atmospheric climatological 

data. In this model, qT. is the temperature at the top level, 

20 meters below the surface. The model is a coarse-grid, 
high=-viscosity model designed to show large scale ocean 


circulation. Complete documentation of the model is given 


in Haney (1974). 
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B. MID-OCBEAN DYNAMICS EXPERIMENT MODEL 


The mid-ocean dynamics experiment, MODE, model is a 
modification of the previously described model. It is 
designed to resolve mesoscale eddies with a horizontal 
grid resolution of 40 km and 5 levels in the vertical. A 
2000 km square section of the general circulation model area 
is represented by the MODE model, located along the western 
boundary and centered at 435N. Circulation is driven and 
controlled by the wind stress and thermal forcing, the same 
as in the previous model, but the eddy viscosity and 
conductivity were reduced to achieve mesoscale eddies. 
Eastern and western boundaries are thermally insulated, 
no-sliv walls, as in the general circulation model. Northern 
and southern boundary conditions are taken as free slip, 
zero normal velocity and a vertical temperature distribution 
which is independent of longitude and time. The ocean 
floor is thermally insulated as in the general circulation 
model. Initial conditions for the model were zero velocity 
and a temperature of 6 at all levels. Twenty years of 
integration was performed i a 100 km X 100 km grid using 
eenoriZzoival viscosity of An = 2x10¢ cm/sec and heat 
diffusivity of Ay = 1x10¢ cm/sec. Ten years of fine grid 
6 


calculations were then made using A. = 2x10 cm“/sec and 


Ay = 5x10° 


om@/sec. 
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TiI. DERIVATION OF AVAILABLE POTENTIAL BNERGY EQUATIONS 





A. DOMAIN AVERAGED AVAILABLE POTENTIAL ENERGY 


The vertically averaged available potential energy 
contained within an idealized ocean basSin may be calculated 
as follows. Assume the basin to be flat bottomed at a depth 
H. The total potential energy per unit volume in a column 


in the basin is given as 


- m 
psd f egzdz , 
Hu 


where Z increases upward from the bottom where Z-=- eH. All 
variations in H are neglected. Neglecting also the dependence 
of density on pressure and salinity, the density may be 
written as a linear function of temperature, equation (6). 


The total potential energy per unit volume in the column 


is then given as 270 T(o) 4 
p= CONSTANT - s«@ | “p A(z) a 2 gr]. 
K 27-H (+4) a 


By defining Z2(T) = 0 if T 2 T(O) and Z(T) = -H if T & T(-H), 
the vertical average total potential energy per unit volume 


in the column may be written as 


Tmax 
P : CONSTANT - g«C, T(-H) 4" + auc. f Zz dT, 
\s 2 H Tie 2 


where Tmin and Tmax represent the coldest and warmest 


temperatures found in the basin respectively. The vertical 
average total potential energy averaged over the horizontal 


domain, denoted by ( ), is given by 


17 





7 - Tax 
P= ConstANT ~ $XG T(CH)H* | anG. f Pe ald) 
A 2 Pe ae ee 

H a 


Traw 
where the bar over Z° denotes an average over an isothermal 
surface. The average potential energy is a minimum when ze 
is a minimum, ie. Z is equal to % everywhere. Writing the 
height of a temperature surface aS a mean height, 2, and 
a departure from the mean, Z', then the domain average total 
potential energy iS a minimum when Z' is zero everywhere. In 


me 
this case, a can be replaced by Z. The minimum domain 


averaged total potential energy can then be written as 


Tmax 
a, = ae 2 => 
Pay = CONSTANT = GAC TOW OH t TX Co Z drt. 
H a H < a 
Moin 


If the difference ‘between the average bottom temperature, 
T(-H), and the average bottom temperature in the state of 
minimum domain averaged total potential energy, T(-H)min, is 
assumed small and thus can be neglected, then the domain 


average available potential energy can be written as 


_ Tmax _ 
A = See f z-# ) dt 
H J 2 } 
Which can be written as 
Cy) 
A a Gf 173 di. (7) 
‘a! 2 
Tin 


Thus, the available potential energy is proportional to the 


variance of height of an isothermal surface. 





B. APPROXIMATION TO THE DOMAIN AVERAGE AVAILABLE POTENTIAL 
ENERGY 


The variance of height over an isothermal surface, gre, 
can be related to the variance of temperature over a horizontal 
Surface, assuming the ocean is stably stratified. When the 
perturbation, &G', from the mean height of an isothermal 
surface is greater than zero, the temperature of this region 
is less than the mean temperature by an amount T'. Similarly, 
when Z' is less than zero, T' is greater than zero. Therefore, 


changes due to Z' are partially compensated by changes in T!. 


The height of a temperature surface may be approximated by 


Z(t) = £(T) =2(1-1), 

which shows the height perturbation can be approximated by 

é t = 

ees “Tez 
OT 

When this is substituted into the expression for the domain 
average available potential energy, equation (7), the 
approximation to the domain average available potential 
energy is then 


Twa 


=\4 
A = ate / He v'dz) dT , 
i 2 oT 
Traiw 
Which can be written as 
, ° ya 
A = 9«C. f {eZ . (8) 
H aT 
oz 
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C. THE BALANCE OF AVAILABLE POTENTIAL ENERGY 


pince available potential energy is a very important 
source for kinetic energy, a study of the energetics of an 
ocean circulation should include a study of production and 
dissipation of available potential energy. Under adiabatic 
conditions an increase in kinetic energy is compensated by 
a decrease of available potential energy. When both energy 
forms increase, then formation of available potential energy 
must be occurring through diabatic processes. A study of 
these generation and dissipation processes and their 
importance may be made through use of the equation of the 
balance of available potential energy. The equation is 
derived as follows, based on the approximate form (8) of 
available potential energy. 

From equation (8), the time rate of change of available 
potential energy is 

aR. nel (2 BT he — ( he [att 2 (Bi). (9) 
oT o(t aT\* ot \oz 
“ (5) & (3) ° 

By use of the thermodynamic equation (5), equation (9) was 


rewritten to show the various terms responsible for the 





change of available potential energy. The term S(t") Sail 


aT" rot , 10 
o(ev} = 1a (10) 


The term 2 (2) can be written as 
at\az 
OT 
° 1 
<i (11) 


be written as 








The local rate of change of temperature and its average were 
obtained from the thermodynamic equation 


> poet ) 2 
of _V- (vt) ~ 2) A AST 4 So SON CED 


The horizontal average of (12) gave 


oT - _ dw) ,KS&T we §(T) | i) 
BY By: 2" 


where the terms -\7-(VT) and Aa7 t were zero Since advection 
and diffusion of heat vanish at the boundaries of the domain. 


Equation (11) was then written using (13) as 











> /dT\). IT dw YT , $M 
s($z)* sel- 3 i ag Js iy 
Which can be rewritten as 
$(2) = 2 |- me 24 i>) 
otloz/° ZL OF ec] ? 
where 
Q _ KT i). 
ec ozL 


Equation (12) was written in terms of an average 
Equation (13) 


temperature and the departure from the average. 
was then subtracted from the result giving 

aT M07) _ Lor), HT), Fr) 977) ee 
= St DE + 3 (USS) ' 3 c (16) 


Equation (16) was then substituted into the right hand side 


oT -_y.Wi)-vqt)- 


of equation (10) to give 


an 





3 (us): Ff ret) - Bevis of vcore] Tle St\ fo ot) «Sey 
26.7 dz J aoe rv(edn) (Ay +56], (17) 


since 


+ OA(wt) - o, 
OZ 


' The first term on the right side of (17) can be written as 


Tl-v-(vt)- Ww) |; -—TWwYVT - 1 TVV ~Twot oT TOW 
Oz oz OZ 


The term -TV:V T is equal to zero because VY T is zero. By 


use of the continuity equation (4), the equation reduces to 


“1. 
4 
< 
-{} 

5 
CV 
£ 

Se 

$3 
~; 
€ 
7d 
+} 
S 


Also, the second term on the right side of equation (17) 


can be written as 


Tv 7 ~ teat) — 2-(¥2") , (19) 
oz | 2 = fee \ 2 , 


The last term of equation (17) can be written using the 


definition of Q as 


FT (Aas) + CB) s S07] : rg (20) 


6 
Substitution of equations (18), (19), (20) into equation (17) 
and using the resulting equation along with equation (15) 
results in the following form for equation (9), the balance 


of available potential energy 
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oz 


( 
Lay - 3a) e a, 


The balance of available potential energy is due to 


(21) 


/ 


five terms on the right side of equation (21). The first 
term is the transformation between available potential energy 
and kinetic energy due to the correlation between temperature 
and vertical velocity. The second term is a triple 
correlation term due to the convergence of the vertical 
transport of temperature perturbations. The third term is 
the generation of available potential energy due to a spatial 
correlation between heating and perturbation temperature. 

The last two terms represent the change of available potential 
energy due to a change in the mean static stability. Lorenz 
(1955) assumed the static stability was constant in time and 
also neglected the triple correlation term leaving only the 
first and third terms on the right side of equation (21). 

In the quasi geostrophic theory, the first and third terms 
are the only terms to appear in the balance of available 
potential energy. For the general circulation model, all 
five terms were studied and for the mid-ocean dynamics 


experiment model the first three terms were studied. 
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IV. FINITe DIFFERENCE EQUATIONS 


The numerical program developed for calculations of 
available potential energy and the balance of available 
potential energy was designed for the general circulation 
model as developed by ees (1974). Therefore, the space 
finite differencing scheme and the boundary conditions 
were the same as those of the general circulation model, 
Haney (1974). A staggered grid system was used for placement 
of the variables. Figure 1 shows the variable placement 
in the horizontal and figure 2 shows the variable placement 
in the vertical. Longitude positions are designated by 
the subscript i, latitude positions by the subscrint j and 
vertical positions by k, increasing downward. 

Horizontal averaging, denoted by q for some variable 
one was done in the following manner. Since the lateral 
boundaries were located at integer i,j points, the mean 
value of a term to be averaged was calculated at the its, jt+d 
positions based on the terms value at, and the areas of, the 
eelaspornvbs at the four surrounding i,j positions. These 
mean values were then summed with equal weights to determine 


the domain horizontal average. 
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A. NUMERICAL CALCULATION OF AVATLABLE POTENTIAL ENERGY BY 
HEIGHT PERTURBATION 


To calculate the available potential energy by equation (7) 
it was necessary to calculate the temperature at the boundaries 


of the layers represented by the subscript k. T represents 


ijk 
the vertical mean temperature for a layer of thickness Da, . 


where 


The temperature and height of the water column at the 


boundaries of each k level were calculated such that Ta 5k 


was maintained as a vertical mean. The actual temperatures 
95 sed? were calculated using the following equations 
y = ~ Teg, - ©2532 ° OZ 
O23 % = a: Tes, ' O4, Oija Like 
DzZy, 


Oish-y = 2 tere * Tin | Ror 2a t KM 


Oriwmeye 2: Tiinm » Okun” Oc inn “DZ km-¥a , 
DZ KM+y, 


where KM represents the total number of levels in the model. 
The horizontal average of the vertical mean temperature, T. ) 


was calculated as explained previously. The depth of oe 


was determined using 95 iked and 44 by use of the equation 


Z(Tk)ijh = Ly-y, pe Sette * Tire (Za (2 ~E rh) 


Ors” Seren 


where 2(T, ) the depth at which the temperature 6 is 


ijk ~ 
equal to the horizontally averaged vertical mean temperature 
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in the layer k. The horizontal average of the depth of 
each of the mean temperature surfaces is A(T)» and 


represents the average height of the isotherms between 8. 
and ®,.1- The departure from the horizontal average was 
determined as the difference between the temperature surface 
ZT) 5 sk and its horizontal average A(T) The resulting 
finite difference analog of (7) for available potential 


energy was of the form 


a _ -_ 
f= See ST Bie (Ba, -Sen)) 
ez 4 


B. NUMERICAL CALCULATION OF AVAILABLE POTENTIAL ENERGY BY 

TEMPERATURE PERTURBATION 

Equation (8) shows the analytical forms of the equation 
for the approximation of average available potential energy 
using the temperature perturbation, T', instead of the height 
perturbation, Z'. 

The horizontal average temperature, T, was calculated at 
each k level. By subtracting this value from the temperature 
at each grid position, a value for the temperature perturbation, 
T', was obtained at each grid position. Using these 
perturbation temperatures, the term ype was calculated at 
each grid position and averaged to the half-integer points 
Greryegts, krt). 

The term for the vertical derivative of the horizontal 


average temperature was calculated at the i,j integer points 
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as 





oT AT Tk — Ty 
+ /. kry = a a 5 
and averaged to the half-integer points. The approximation 
to the mean available potential energy was then obtained by 


combining these terms and summing over all levels as indicated 


KM 
eee, 


re 


by 





Coe 
Ma, : 


a 
£4 heey Zirh (22) 
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Kth, 


C. NUMERICAL CALCULATION OF THE BALANCE OF AVAILABLE 
POTENTIAL ENERGY 
The analytical form of the rate of change of available 
potential energy is shown by equation (9) as 


a c | faa3 iv \oe -( ee 2 (3) - (9) 
vs (33) ye (3) Ot \oz 
OZ 


pre ana Yea in (9), were calculated as in 


The factors i212 
the numerator and denominator of (22). The remaining parts 
of equation (9) were calculated by use of the right side 
of the thermodynamic equation as indicated by equations (10) 
and (11). Equation (10) was written to determine 2.(47") 
ey ecatlenilating Tot at integer values of i,j,k and then 
averaging to obtain its value at half-integer points. 
Equation (11) was calculated from 

> (2%) : OT ci k = OT Ej he 
oz Ob Aiken a se ae 
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and then averaged to the half-integer grid points. The 
resulting finite difference form for equation (9), the 


balance of domain average available potential energy, is 


then 
- Ke- | : 
dA - 3&5 Fe ° 2°") ~ (%t"} 2 fot | 
——~ eee =7 ease | Ge enn ORF mee | ee ee (23) 
at tH < (2 | >* 2 (22) OZ \ ot 4 
i h a oz Reh we 


D. NUMERICAL CALCULATIONS OF THE TERMS FOR THE BALANCE 
OF AVAILABLE POTENTIAL ENERGY 
As indicated by equation (21), the balance of mean 
Pees labile potential energy contains five terms which are 
Beeained by use of the thermodynamic equation cS). 
Development of equation (21) showed that the five terms 


could be written as below 

















“Twat : T [9 (uF) = ACR] (18) 
8) Ew 
rs (ast) + (Kz) ray] (20) 
af. 2G ): a(- an) (24) 
2 2.) can o\ ee + —_) (25) 








The terms in the brackets on the right side of the above 

five equations were calculated using the identical finite 
difference method as used in the general circulation model 
to calculate the terms of the thermodynamic equation. The 


details are in Haney (1974). 
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To make an accurate comparison of the two expressions 
for domain average available potential energy, equations (7) 
and (8), and to determine the importance of the terms 
occurring in the balance of available potential energy, 
equation (21), the general circulation model, as described 
in chapter II, was run for large-scale motion. A description 


of the conditions and the results obtained are Given below. 


A. GENERAL CIRCULATION MODEL OF THE OCHAN 


To establish quasi-steady state conditions, the model, 
Haney (1974), was run in two phases. For the first phase, 
the linearized form of the vertical shear current was used. 
The vertical shear current equations were obtained from the 
primitive equation form of the horizontal equations of motion, 
(1) and (2), by subtracting the vertical mean of equation (1) 
from equation (1) and by subtracting the vertical mean of 
equation (2) from equation (2). The vertical mean current is 
assumed constant in time. The temperature is determined from 
the complete thermal equation ‘5). The first phase integration 
was for approximately 100 years. For the second phase, the 
complete nonlinear equations were used. The second phase 
integration was for 25 years. A more complete description of 
phase 1 and phase 2 may be found in Haney (1974). Twenty-four 


sets of model history data was taken at this point at intervals 


De 








of 25 days, representing a total of 600 days. The results 
showed the motion to be not strictly steady-state. As shown 
by Haney (1974), slight warming was occurring near the bottom 
after the integration period indicating a continuing downward 
flux of heat. Figure 3 shows the temperature at 20 meters 
below the surface, averaged over the last year of integration 
for phase 2. Figure 4 shows the vertical velocity at 60 meters 
below the surface averaged over the last year of the second 
phase. Comparing figures 3 and 4 shows two large gyres 
centered on each side of the equator. In these gyres, warm 
water is being forced down. Near the equator and along the 
northern and southern boundaries, cold water is being forced 
upward. Thus, available potential energy is being generated 


vy these motions. ae 
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Ficure 3. Temperature (°c) at 20 m, averaged over 
tne last year of the second phase. 
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1. Results of Available Potential Bnergy Calculations 


Comparison of calculations of available potential 
energy as determined by equation (8) and the approximate 
form (9) Showed that the available potential energy as 
calculated by the height perturbation to be slightly larger, 
(1383 ergs/cm) , than the available potential energy as 
calculated by the temperature perturbation, (1232 ergs/cm”), 
figure 5. This discrepancy is perhaps due to Tmax and Tmin 
being not uniquely defined in Z coordinates. Fluctuations 
in available potential energy-due to a minor change in the 


model near the end of the 600 day period of record were 
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Figure 5. APE by #Z'“(e) and by 31'“(e) for the 
General Circulation Model. 
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Simulated by both forms of available potential energy 
calculations. Figure 5 indicates equation (8) calculated 
available potential energy to a sufficient degree of accuracy 
to be considered a good approximation equation for available 
potential energy in this model. 
e. Results of Calculations of the Balance om 

Available Potential Rnicisas 

From equation (21), the balance of available potential 
energy is represented i une sum of two integrals, 


Di = heater 


Ot 








where Sea — 
parta = 3% - we to (v}« AL ake 
ae (23 az \2 (> oe 
" 1. oe (eke | oe (26) 
part = 3AxC 4T't 2 [L3 a(wr) , 9. dz. 
pa Tir ary 32 32 "ec 
& 4 5. 


According to their definitions, Q! and Q contain the convective 
adjustment of temperature, $ (t)!, and 8 (7) respectively. 
In the calculations below, these convective adjustment terms 
were omitted due to the difficulty of recovering this 
irreversible process from the history records. 

Part 1 represents changes in available potential 
energy resulting from changes in the variance of the temperature 
perturbation on a horizontal surface. Heating of warmer 
regions and cooling of colder regions at the same depth 
increases the horizontal temperature variance and increases 


available potential energy. Part 2 represents changes in 
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available potential energy due to changes in mean static 
stability. Heating of upper levels and cooling of lower 
levels increases the static stability and decreases available 
potential energy. The magnitude of these integrals are 

shown by figure 6. Variations in these two parts after 

day 400 of the history record were caused by the fact that 
the convective adjustment, SCE was changed to double 
precision accuracy in the model. 

The processes represented by part 1 are Serving to 
decrease the total available potential energy in this model, 
indicating the model, apart from d(T) ', is working to 
reduce temperature perturbations. The magnitude of the 
three terms comprising part 1, equation (26), are shown 
Dieitegunes,, where all three tenms are of the same order of 
magnitude. The first term of part 1 shows the transformation 
of energy between available potential energy and kinetic 
energy due to a correlation between the temperature 
perturbation and vertical velocity. Results indicate that 
in the domain average, energy is being transformed from 
kinetic energy to available potential energy. This is 
consistent with the results of Haney (1974), as discussed 
above, where a transformation of kinetic energy of vertical 
shear motion to potential energy was occurring. Term 2 
represents a triple correlation due to convergence of the 
vertical flux of temperature Puce eens. Being a triple 
correlation of perturbation quantities, Lorenz (1955) 


assumed this term to be negligible in the atmosphere. In 
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this general circulation model, term 2 is smaller than all 

of the other terms by only a factor of 2 to 3’'and thus it cannot 
be neglected. This term arises because the static stability, 
OT/OZ, is not constant in Z. If it were, the term would 
vanish. The term would also vanish if the triple correlation 
were small. One reason it was not small is because of the 
forced nature of the vertical velocity, w, and temperature 
fields near the boundaries. With heat flux through the 
boundaries, neither the vertical velocity nor the temperature 
gradient in the horizontal would have been as large as 
occurred. The third term of part 1 shows generation of 
available potential energy due to a spatial correlation 
between heating and onenneaion temperature. This term 

was calculated without the contribution from Sa). 

present in equation (20). The negative result of term 3 
indicates a loss of available potential energy due to the 
non=convective heating. The magnitude of this term indicates 
the importance of the thermal diffusion terms in this 

oceanic general circulation model. 

The processes represented by the terms of part 2, 
figure 8, were of the same order of WeeRanee. and tend to 
increase the available potential energy in the model. The 
Mirstesern Of pact ve, term 4) as due Co the vertical flux 
of mean temperature and the second term, term 5, is due to 
mean heating. As in term 4, the sonveetive adjustment 
tern, a(t), was absent from the mean heating term. 9 (T) 


tends to increase the mean static stability, OT/OZ, and to 
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decrease available potential energy. It's inclusion would 
allow terms 4 and 5 to sum to a negligible value. Lorenz 
(1955) assumed part 2 was zero for steady-state atmospheric 
motion. Under conditions of seasonal radiation heating, 
the vertical flux of mean temperature, term 4, may not be 
in phase with the mean heating that occurs, term 5, and 
thus would not cancel. The general circulation: model was 
not designed to show seasonal radiation heating. 

The energy flow diagram depicting terms 1 - 5 in 
equation (26) is shown in figure 9. Terms 1,4 and 5 were 
serving to increase available potential energy by an amount 
of 1.705x107? ergs/cm’sec. Terms 2 and 3 were serving to | 
decrease available potential energy by an amount of .908x1077 


d 


ergs/cm’sec. The available potential energy during the 
record period was observed to be decreasing, figure 5, with 
an average value of 1232. ergs/cm”. The discrepancy between 
the terms serving to increase available potential energy, 

Maree | = part 2)0, and the observed decrease in available 
potential energy, figure 5, was determined to be due to not 
including the é.(f) term in the balance of available potential 
energy calculations. The value of the term required to 
balance the energy flow diagram, as in steady=state conditions, 
is indicated by term 6 and shown by the dashed line in 

figure 9 to be approximately 7972x107 ergs/cm’sec. These 
results indicated convective adjustment was an important 


process in the model, as was expected since the model 


contained a large region of near isothermal water at high 


42 





Tete \.oeo) 









Terms ( s/c) 


APE 
(1232) 


“Perm (.610) 





Term2 (.233) 


vw Term6 (.797) 
Terml (.623) 


UNITS: 
APE.......2...ergs/em- 
TERMS i-6...10 "ergs/cm”’sec 


Figure 9. kLnergy Flow Diagram for the 
SeuecroieGirewarion Hodel. 


4 





latitudes where convective adjustment would frequently 
occur. .The convective overturning plays no role in the 
kinetic energy budget and therefore the arrow for term 6 


Shows this effect to be an external loss to the systen. 


B. MID-OCEAN DYNAMICS EXPERIMENT MODEL 


Figure 10 below shows the instantaneous streamfunction 
for day 10080 during the seventh year of fine grid 
calculations. Twenty~four sets of data for the experiment 


were taxen from day 10050 to 10110 at 23 day intervals. 





Figure 10. Instantaneous streamfunction for 
day 10080 of the Hid-Ocean 
Dynamics Experiment Model. 





To determine the available potential energy as a result 
of transient motion and standing motion, equation (8) was 
written with temperature in terms of a time average 
temperature, ([T], for standing motion, and a departure from 
the time average, T, for transient motion. Since available 
potential energy is determined from the temperature 
perturbation, the horizontal average of the time averaged 
temperature was determined. Then the perturbation temperature 
of the standing motion, [T]' and transient motion, rT, was 


determined as below 


PiedilesT" 
oe eo oe ee (27) 


where {T] is the horizontal averaged standing temperature 
caer” is the horizontal averaged transient temperature. 
The term [T] '! is the standing eddy temperature and T. is 
the transient eddy temperature. The eddy temperature is 
then the sum of the last two terms of equation (27). 
Substitution of the eddy temperatures into equation (8) 
resulted in the following form for domain Beret available 


potential energy, 
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The time averaged available potential energy was written 


as 
° , —=any TF 
PE 2 xt ff {Tl nie J} a 
eae + +t Ge 
su 


* 1 : 
where the term \s(ry'2 | vanishes. In this model, the 


C25) 


static stability factorLla T/aZ| was assumed constant in 
time and written as ® . Equation (29) shows time average 
available potential energy as the sum of a part due to the 
standing field and a part due to the transient motion. 

For the study of the balance of available potential 
energy in the mid-ocean dynamics experiment model, only 
the first three terms of equation (21) were considered 
important for further study. Because lateral boundary 
conditions allowed heat to flow in or out of the model, it 
was reasoned that convective adjustment would not play an 
important role in the model. With little convective 
adjustment and steady forcing, it was therefore assumed 
that the long-time average of the sum of terms 4 and 5, 
part 2 of equation (26), would be negligible. Calculations 
of the balance of total available potential energy, below, 
showed the terms of part 2 summed to approximately zero 
for each set of data as well as in a time average so the 


assumption was justified. The study therefore concentrated 





on transformations of energy due to coe: From the first 
three terms of equation (21), the balance of available 


potential energy was approximated as 
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Equation (30) was written in terms of a part due to transient 
¥ 
motion, q , and a part due to standing motion, (qj; for a 


quantity q and then time averaged as 


[24] = Sete ((- fo) [1 Lt)" ~Twrr® | at i eld Cr]',.0w* twtr] 
a} oH 4, ; % vba 2 
Hon _ 1b. 2a. 2b. 


<a, - mee) oer cen 
2 v 8 
Ow ed. 4a. Dior 
where the quantities containing only one term due to transient 
motion were zero in the long-time average. The values of 
terms ec and ed were not calculated in this study but their 
Sum could be inferred from terms 2, 2a and 2b. Since heat 
diffusivity was allowed across the boundaries, the definition 
of Vc in (25) was changed to the following 
od nathan signe omnia ———— 
2, % 
@ 2 (ATT) . k ot) + Se(T) 
Ee yA 


Since OS da waS not equal to zero in the MODE model. 


1. Results of Available Potential Energy Calculations 


Figure 11 shows the total available potential energy 


as determined from equation (28) to be slowly increasing 


during the history record. The first term on the right side 
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Figure 11. APE vs. Time for the MODE Model 
(Total Motion) 
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of equation (28), due to standing motion was constant in time 
and had a value equal to 116.65 ergs/cm’, shown in figure 11. 
This was an order of magnitude smaller than the available 
potential energy of the general circulation model and it 
points out the importance of the horizontal diffusion of 

heat across the northern and southern boundaries in preventing 
large horizontal gradients of temperature. The.third term 

of equation (28), due to transient motion, is shown in 

figure 12. This term apvpears to have a period of approximately 
100 days, but is too small in amplitude to significantly 
affect the total available potential energy. Since the left 
Side of equation (28) was known from figure 11, and two 

terms on the right side of equation (28) were known, the 
middle term on the right side of equation (28) was found to 

be primarily responsible for the slow increase of available 
potential energy seen in figure 11. The time average of 

the total available potential energy, corresponding to the 
left side of equation (29) was 117.3 ergs/em’. The first 

term on the right side of equation (29), the available 
potential energy due to standing motion, is the same as the 
first term on the right side of equation (28), equal to 

116.65 ergs/cm?. Thus the time average available potential 
energy due to transient motion, the last term of equation (29), 
had a value of .66 ergs/cm? and was only a small part of the 
time average total available potential energy. Being so 
small, the available potential energy due to transient 


motion would probably not be noticed on a temperature map 
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Figure 12. APE by #f'° vs. Time for the MODE Model 
(Transient Motion) 
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from the model. If the variations in kinetic energy were of 
this order of magnitude, it could be important in the kinetic 


energy balance which would be subject to future investigation. 


2. Results of the Balance of Available Potential 

Knergy Calculations 

By the use of flow diagrams, the magnitude of the 
terms responsible for the balance of available potential 
energy can be easily seen. Figure 15 shows the balance of 
the time average total available potential energy, 
corresponding to equation (30). This figure indicates energy 
was being transformed from kinetic energy to available 
potential energy by term i of equation (30) at a rate faster 
than it was being removed by terms 2 and 3. The average 
rate of increase of total available potential energy, as 


3 


shown in the diagram, was approximately 046x107? ergs/cmsec. 
From the observed change in available potential energy, 
figure 11, there was an average rate of available potential 
energy change during the 60 days of .05x107? ergs/cm SCC. 
Figure 14 shows the value of the energy transformations 
due to standing motion, corresponding to terms la, ea, 43a of 
equation (31), and figure 15 shows the value of the energy 
transformations due to transient motion, corresponding to 
terms 1b, 2b, 3b of equation (31). The average rate of 
increase of available potential energy by the standing field, 
as shown in the diagram, was approximately 075x107? 


3 


ergs/cm” sec. The average rate of decrease by the transient 
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field, as shown in the diagram, was 035x107? ergs/cm” sec. 
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Instead of such changes in the available potential energy 
fields, this perhaps simply indicates that a conversion of 
available potential energy from the standing to the transient 
fields was occurring. The conversion term for available 
potential energy between the standing and transient fields 
was not calculated because it also contained terms of 
Standing and transient motion as does terms ec and ed. The 
conversion term must have had a value between .055 to .075 
x107? ergs/cm? SCC. Summing of the net increases by the 
standing field and the transient field gave a net increase 
for the total of .040x107? ergs/cm? sec, which compares 
with the value of 046x107? ergs/cm? sec seen in figure 13. 
The sum of terms 3a, figure 14, and 3b, figure 15, was equal 
to term 5 of figure 15 as it should be, since in the time 
‘average term 3 is composed of only term 45a and 3b. Similarly, 
term 2a, figure 14, and eb, figure 15, sum to equal term 2 
of figure 13. However, from equation (31) term 2 is composed 
of four terms in the time average. Therefore, terms ec and | 
2a must be negligible or sum to a negligible value. 

The sum of 1a and Ib is less than term 1 by an 
amount of 006x107? ergs/cm? sec. some round-off errors 
may have occurred to cause the small discrepancy. Comparing 
figure 143 and figure 14, the total energy being transformed 
from kinetic energy to available potential energy, term 1, 
was almost totally by the standing field, term la. Since 


term 1b, figure 15, was small and positive, the source of 


kinetic energy for the transient motion was not available 


DE, 





potential energy during the 60 day period of the MODE model. 
Figure 16 shows the total transformation from kinetic energy 
to available potential energy, term 1, was undergoing large 
fluctuations. At some times the transient eddies did draw 
energy from available potential energy, since ib, figure 17, 
was at times negative, but the energy flow in either direction 
was small. These results may indicate that the. transient 
eddies were not baroclinic for the 60 day period. The 
amplitude of the eddies of the MODE model were known to be 
large. With large amplitude eddies, the dissipative forces 
were expected to be large, indicating the transformation of 
energy from the source of energy for the kinetic energy of 
the eddies was large. Therefore the primary source of energy 
for the transient eddies was not available potential energy 
during the 60 day history record. 

It was suspected that errors may have occurred due 
to the 23 day rate of sampling the data. Figure 16 shows the 
three terms of equation (30) as a function of time. Term 1 
and 2 fluctuate to a greater extent than term 3. Since 
both term i and 2 are correlated with vertical velocity, w, 
the possibility of aliasing is further supported. Fluctuations 
in term 1b, figure 17, and 2b, figure 18, also lends support 
to this possibility due to their fluctuations. Neither 
term 3 of figure 16 nor term 4b of figure 19 showed a large 
degree of fluctuation and neither of these terms were 
correlated with vertical velocity. Therefore, the above 


energetics were recalculated using a 120 day history record 


56 





0.4 


eT \ / [ / | 


eres 
ae 0.1 


0.0 
-0.1 
% RH HH RM # % HH 
eR HK HH ae =: z 
WH xX 
-0,2 7 a eG 
a = \ 
Lf J oe i " a 
8 
-03 \ e \ / > \/ 
s a & . 
8&8 | | 
a CG 
-04 
1.0 30.0 60.0 


DAY of RECORD 


Figunes|G. Contribution by Term1x107(e), 
MepmexlO°(e), Term3xl07(x), to the 
Balance of APE, Eq. (30). 


De 








1.0 30:9 60,0 
DAY of RECORD 


Figure 17. Contribution by Termib x10? 
to the Balance of APE. 
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. Figure 18. Contribution by Term2b x104 
to the Balance of APE. 
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Figure 19. Contribution by Term3b x10 
to the Balance of APE, 
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which began on day 10050 and ended on day 10170. No essential 
changes occurred in any of the terms except term 1b, which 

was still small, but negative for the 120 day average. Since 
a 120 day average would probable be more accurate than a 60 


day average, the eddies were felt to be baroclinic. 
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VI. CONCLUSION 


The results obtained from this study indicated that 
both the general circulation model and the MODE model were 
mechanically driven by the wind stress and not thermally 
driven to a great extent. The below diagram shows the 
transformation of energy from its generation by the wind 
stress to its dissipation through transient eddies as 


indicated by the results from the MODE model. 





The atmospheric wind circulation acts on the ocean surface 
creating large scale standing motion. Some of the kinetic 
energy of tne standing motion is transformed roms hot beom mets. 
available. potential energy of the ocean. Approximately 
half of the input of energy to standing available potential 


energy is lost through non-convective heating by standing 
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motion. About one-fourth of standing available potential 
energy is transformed to transient available potential energy, 
where a large portion is lost through transient non-convective 
heating. A small portion, in the transient available 
potential energy balance, is transformed to kinetic energy 

of transient motion where it is lost through dissipation. 
Whether the transient available potential energy is the 
primary source for the transient kinetic energy must be 
determined by a study of the transient kinetic energy 

balance. 

In the atmosphere, forcing is by the heating and the 
energy transformations are from available potential energy 
of standing motion to transient available potential energy, 
then to transient kinetic energy. 

The main conclusion is that baroclinic instability is 
@ possibility for maintaining the MODE eddies. The general 
circulation model showed energy transformations from eddy 
kinetic energy to eddy available potential energy for large- 
scale quasi-steady state motion. The MODE model showed the 
energy transformations to be opposite, from eddy available 
potential energy to eddy kinetic energy, for transient 
motion. Therefore, during periods of baroclinic instability, 
eddy available potential energy is aiding in the maintenance 


of the eddies. 
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